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ABSTRACT 


The performance characteristic of a digital correlation 
radiometer signal processing back end (DBE) is analyzed 
using a simulator. The particular design studied here 
corresponds to the airborne Hurricane Imaging radiometer 
which was jointly developed by the NASA Marshall Space 
Flight Center, University of Michigan, University of Central 
Florida and NOAA. Laboratory and flight test data is found 
to be in accord with the simulation results. Overall design 
seems to be optimum for the typical input signal dynamic 
range. It was found that the performance of the digital 
kurtosis could be improved by lowering the DBE input 
power level. An unusual scaling between digital correlation 
channels observed in the instrument data is confirmed to be 
a DBE characteristic. 


Index Terms— digital radiometer, HIRAD, digital 
correlation, kurtosis, correlator. 


1. INTRODUCTION 


In modern microwave radiometer systems used for 
interferometric, polarimetric and/or spectrometric 
measurements of natural thermal radiation, digital signal 
processing offers many practical advantages [1]. The 
interferometric measurement requires signals from many 
different antennas to be correlated whereas signal from the 
same antenna at two different polarizations is correlated for 
the polarimetric measurements. For spectrometric 
measurements one option is to correlate the antenna signal 
with a time-delayed version of itself to map out the auto- 
correlation function. Then by Fourier transformation the 
power spectrum can be obtained. The other method is to 
divide the radiometer passband into multiple subbands using 
polyphaser filter bank for spectral measurements. Polyphase 
filter banks can also be used for radio-frequency interference 
(RFI) mitigation [2]. In many modern radiometric systems 
combinations of these techniques are implemented. 
Furthermore digital correlation calls for analytical signals 
(Inphase, I & Quadrature, Q signals) to be derived from the 
received antenna signal. This can also be implemented 
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Figure 1: HIRAD system block diagram with DBE 


digitally using Hilbert Transform (HT). Due to the advances 
in the development of modern high speed FPGAs (Field 
Programmable Gate Array) different signal processing 
schemes can be implemented and investigated using the 
same hardware. This offers to a great amount of flexibility. 

Digitization involves sampling of the voltage signal 
at periodic interval and quantizing the sampled values so 
that each can be represented by a finite number of bits. The 
quantization effect on digital correlation has been studied 
extensively [3-13] but these studies does not include end-to- 
end effect of finite bit truncation along the Digital Signal 
Processing (DSP) chain, if further processing of the 
quantized signal involved before computation of the digital 
correlation. For analytical signals, there may be a gain 
imbalance between I and Q channels because of HT filtering 
and bit truncation [14] which might also affect the 
correlation between quadrature channels. Here we present a 
study of a radiometer digital back end (DBE) which digitizes 
a band limited Intermediate Frequency (IF) Gaussian noise 
signal, performs polyphase filtering, I/Q demodulation and 
finally computes cross-correlation and statistical moments. 
The particular design studied here corresponds to the DBE 
of the airborne Hurricane Imaging Radiometer (HIRAD) 
jointly developed by NASA Marshall Space Flight Center, 
University of Michigan, University of Central Florida and 
NOAA [15]. 


2. HIRAD DBE 
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Figure 2: DBE non linearity (%) v VS. : eeaeenes (dBm) 
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Figure 3: Gig deviation (%) from 1.296 vs. input power 


HIRAD instrument uses interferometric aperture synthesis 
technique to create high resolution cross-track brightness 
temperature image at 4.0, 5.0, 6.0 & 6.6 GHz using multiple 
fan beam antenna arranged in a thinned array configuration 
[16]. Each antenna is connected to a dedicated receiver 
which sequentially translates four RF passbands to a 
common 75 MHz band centered at 187.5 MHz by switching 
local oscillator (LO) frequency. This Intermediate 
Frequency (IF) signal is digitized by the DBE for further 
processing. The entire DBE is implemented in FPGA and its 
configuration is identical for all four RF frequencies. A 
schematic of a pair of antenna channels in HIRAD system is 
shown in Fig. | with the DBE inside the dashed box. The 
Analog to Digital Converter (ADC) samples the 150-225 
MHz IF signal at Nyquist Rate (150 MSa/s) with 12 bit 
resolution. However, lower 5 bits are discarded and only 7 
bits are transferred to the subband filter stage. The subband 
filtering is done by 16 parallel polyphase bandpass FIR 
filters spanning the digital baseband of 0-75 MHz. The filter 
coefficients are represented using 9 bits and the output is 22 
bit wide. A copy of each subband filter output is fed into a 
HT FIR filter (Q channel) and another is just phase shifted 
to match the phase shift due to HT filter (I channel). HT 
filter coefficients are 9 bit and the filter output is 35 bit 
wide. Therefore I channel is 22 bit wide and Q channel is 35 
bit wide at the output of I/Q demodulation block. Only 7 
most significant bits (MSB) are selected from both I and Q 
channels for further processing where cross-correlation (In, 
QmQn, LnQn & Qml, : subscript m,n will be dropped here- 
after) and moments (1“, 2"*, 3™ & 4'") are computed. This 
method of I/Q processing introduce a digital gain imbalance, 
Giq [14] which needs to be known for instrument 
calibration. Giq is the amount of scaling required for the Q 
channel to match I channel output. The theoretical Giq value 
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Figure 4: deviation (%) of kurtosis vs. input power 


of 1.2975 is obtained by considering the effect of scaling 
difference due to bit truncation and amplification of the Q 
channel signal due to Hilbert filter convolution. This value is 
very close to the value 1.296 obtained from instrument data. 


3. DBE SIMULATION 


A computer program is developed to simulate the HIRAD 
DBE signal processing including ADC quantization, bit 
truncation and finite bit arithmetic. A pair of discrete time 
partially correlated noise signals are generated in software 
whose power is controlled by changing the variance 
assuming a 50 ohm system. The magnitude of the complex 
cross-correlation is set to | and the phase angle is varied to 
simulate different correlation for the real and imaginary part 
of the analytical signal. Thousand ~9ms integration periods 
are simulated for each power level between -30 dBm to +1 
dBm in one dBm steps. It should be noted that the output 
power of HIRAD receivers were measured in the lab ( @ 30 
deg C) for all possible antenna switching states spanning the 
entire receiver dynamic range from cold load (min) warm 
load + noise diode (max). The dynamic range is between -5 
to -2 dBm which included all 4 possible LO frequency. To 
verify the working of the simulator, the software generated 
noise signals were input (using an arbitrary waveform 
generator) to a spare HIRAD DBE running the same FPGA 
logic as the flight system. The FPGA recorded data agreed 
very well with the simulated data. 


3.1. Moments vs. Input Power 


The most important statistical moment of the signal is the 
second moment which is proportional to the power of the 
signal under measurement. Therefore the digital second 
moments from I and Q channels are a measurement of the 
system noise temperature from which the antenna 
temperature can be determined. The I and Q_ second 
moments from the same antenna are redundant 
measurements which can be combined to_ reduce 
measurement noise. Since there is a Giq scaling between the 
I and Q channels the scaling for the second moments should 
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interfering 5 GHz WiFi signal (top). Same with a 10dB 
attenuation of the IF signal before digitization. 


be Giq’. The 1.296 value of the Giq was derived using 
second moments recorded from HIRAD flight data. In Fig. 
2, the % nonlinearity of second moment counts are plotted 
with respect to the input power for I and Q channel 
separately. Various colors are for 16 different subband 
channels. The HIRAD dynamic range is marked with two 
vertical dashed lines. Note that above -3 dBm % deviation 
starts to deviate from 0. Even at -1 dBm the non-linearity is 
less than 0.1%. The Giq computed from simulated data 
matches perfectly with that observed real flight data except 
subbands #1, #8, #9 and #16 (bands near passband edges 
and center). Fig. 3 shows the % deviation of Giq from the 
value 1.296 for various input power. HIRAD dynamic range 
is in a “sweet spot” where Giq is essentially flat for most of 
the subbands. 

The higher order moments are used along with the 
2™ moment to estimate kurtosis of the input signal. Natural 
thermal emission (i.e. Gaussian noise) has a kurtosis value 
equal to 3. If the radiometric measurement is contaminated 
by manmade radio frequency interference (RFI) then the 
detected signal kurtosis can have any value. Therefore RFI 
which makes signal kurtosis other than 3 can be detected by 
HIRAD DBE. However, quantization affects kurtosis [] and 
its effect in HIRAD DBE is studied using the simulator. The 
simulated % deviation of kurtosis from the ideal value 3 is 
plotted in Fig. 4 for various input power level. Note that the 
kurtosis has strong input power dependence in the HIRAD 
dynamic range. The I channel is affected more that the Q 
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channel because of the relative digital scaling between the 
two. The same phenomenon is also observed in HIRAD 
flight data where the 2"’ moment is found to be correlated 
with the kurtosis. Further, kurtosis showed a variation with 
the antenna switching states. HIRAD receiver switches 
between antenna, warm and cold load. For each of the 
switching state a noise diode (nd) signal is turned on and off 
creating six possible radiometer states. Top panel of Fig. 5 
show the instrument kurtosis data collected in a lab 
environment for various switching states (30 sec average). 
Different colors are various subband channels. As expected 
from the DBE simulation, the kurtosis varies between 
receiver states. The bottom panel of Fig. 5 is plotted from 
data collected few minutes after with a 10 dB attenuator 
placed between the receiver and DBE. All the kurtosis 
values becomes very close to 3 for all receiver states. 
However, for subband #8 and #9 (center of the 5 GHz RF 
passband) the antenna and antenna + nd signal kurtosis 
deviate significantly - which wasn’t readily discernable in 
the data from original HIRAD configuration (Fig. 5, top 
panel). The lowering of the kurtosis value is found to be due 
to interference from the 5 GHz WiFi signal. Therefore, 
kurtosis sensitivity can be improved by lowering DBE input 
signal power as predicted by simulation results in Fig. 4. 


3.1. Cross-correlation vs. Input Power 


An unusual gain imbalance between II (QQ) and IQ (QI) 
counts was noticed during interference pattern measurement 
in the anechoic chamber against a radiated noise diode. 
Because of the Giq scaling between I and Q channel one 
would expect the maximum peak to peak swing of II to be 
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Figure 7: RMS error in digital correlation plotted against 
DBE input power 


Gig? times that of QQ and Giq times that of IQ (or QD). 
However, interference pattern counts from the chamber data 
showed that (top panel Fig. 6) IQ (QI) to be smaller than 
QQ. The relationship between II and QQ still held true. This 
behavior of the DBE wasn’t expected but the simulator 
confirmed that this is a characteristic of the implemented 
DSP. In the Fig. 6 bottom panel the simulator output is 
plotted for various correlation input which show the exact 
same behavior observed during the chamber test. Next the 
normalized digital correlation is computed for various input 
analog correlation at different input power level. The RMS 
difference between analog and normalized digital correlation 
combining II, IQ, QI & QQ is plotted in Fig. 7. In the 
HIRAD dynamic range (between the dashed vertical lines) 
The RMS error is below 0.025 for all subbands. 


4. CONCLUSIONS 


The performance of a radiometer digital back end (DBE) is 
analyzed here. The particular design studied here 
corresponds to the DBE of the airborne Hurricane Imaging 
Radiometer instrument. A computer simulator is developed 
to analyze effect of input power on various DBE output 
products. The 2"’ moment nonlinearity is found to be 
negligible in the HIRAD dynamic range. However, kurtosis 
statistics can be improved by lowering the DBE input power. 
Digital correlation error is found to be minimal in the 
HIRAD dynamic range. An unexpected scaling between 
II(QQ) and IQ(QI) observed during the interference pattern 
measurement is confirmed to be a DBE characteristics. 
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